Multipartite Entanglement Evolution Under Separable Operations 
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We study how multi-partite entanglement evolves under the paradigm of separable operations, 
which include the local operations and classical communication (LOCC) as a special case. We prove 
that the average "decay" of entanglement induced by a separable operation is measure independent 
(among SL-invariant ones) and state independent: the ratio between the average output entangle- 
ment and the initial entanglement is solely a function of the separable operation, regardless of the 
input state and of the SL-invariant entanglement measure being used. We discuss the "disentangling 
power" of a quantum channel and show that it exhibits a similar state invariance as the average 
entanglement decay ratio. Our article significantly extends the bipartite results of 0-0] as well as 
the multi-partite one of [|| , all of the previous work being restricted to one-sided or particular noise 
models. 

PACS numbers: 03.67.Mn, 03.67.Hk 



Introduction. Entanglement is a key ingredient in 
many quantum information protocols, such as factor- 
ing Q, teleportation Q, quantum dense coding Q or 
measurement-based quantum computation Q. In all 
these schemes entanglement plays the role of a resource 
that is "consumed" during the actual implementation of 
the protocol. For example, in measurement-based quan- 
tum computation, a highly multi-partite entangled state 
(the cluster state) is adaptively measured locally until 
the desired unitary evolution is achieved; the cluster state 
model is equivalent to the well known circuit model @, 
but, in contrast to the latter, the role of entanglement as 
a computational resource is now clearly visible. 

Decoherence (To| constitutes a major impediment in 
the physical implementation of quantum information pro- 
tocols. In general, entangled states are not resilient 
against noise, and this may first look like an impassible 
barrier for the construction of a working quantum com- 
puter. Quantum error correcting codes can make quan- 
tum computation possible provided the level of noise is 
below a certain threshold, but the current technology is 
far from achieving it. It is therefore of crucial importance 
to understand how entanglement behaves under decoher- 
ence. 

Previous work [Tj-Q addressed this important prob- 
lem of entanglement evolution (decay) , however the noise 
model was restricted to one-sided channels: i.e. noise 
that acts non-trivially only on one subsystem, without 
affecting the other remaining parts of it. 

In this article we investigate multi-partite entangle- 
ment evolution under a far more general setup: the 
locally-correlated noisy channels: technically separable 
operations, which includes the LOCC as a strict subset 
fill ] . Separable operations (and automatically LOCC) 



constitute a much richer class and encompasses the pre- 
vious investigated one-sided channels as well as the usual 
independent noise model assumed in the theory of quan- 
tum error correction. The study of entanglement evolu- 
tion under such noise model is therefore deepening our 
understanding of the fascinating and not-at-all well un- 
derstood subject of multi- partite entanglement fl2| . 

Our most surprising result is that entanglement evolu- 
tion is independent of the initial state and the entangle- 
ment measure used to quantify its evolution, and depends 
solely on the intrinsic properties of the noise. We provide 
a closed-form formula for how entanglement decays which 
is straightforward to calculate given the Kraus operator- 
sum representation of the process. 

Recent attempts to characterize entanglement evolu- 
tion under locally-correlated noise, not just one-sided, 
were made in [l3[, however the authors restricted their 
study to tripartite qubit entanglement and particular 
form LOCC protocols. Our result is the first one that 
completely characterizes the average entanglement evo- 
lution under the very general class of separable opera- 
tions and may be of interest from an experimental point 
of view. 

Entanglement evolution. Wc consider a multipartite 
system of n qudits, described by a Hilbert space H = 



Hi <g> H 2 



T~L n - and let cL be the dimension of 
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the i-th local Hilbert space. Denote by B(H) the set 
of all bounded operators (e.g., density matrices) act- 
ing on %. We further define the special linear group 
G = SL(di,C) (8 SL(d 2 ,C) (g> • • • ® SL(d n ,C), where 
SL(d, C) is the group of d x d complex matrices of deter- 
minant 1. 

Definition 1. An SL-invariant multi- partite entangle- 
ment measure Q/ } Ei nv (-), is a non-zero function ini- 
tially defined from pure states in % to non-negative real 
numbers then extended to mixed states satisfying: 

i) It is SL-invariant, i.e. Ei nv (g\tp)(i(;\g ii ) = 



2 



E inv (\ip}(tp\), f° r al1 9 £ G and all \ip) E H; 

ii) It is homogeneous of degree 1, i.e. Ei nv (r\tp)(ip\) = 
rE inv (\ijj}(ip\), for all non-negative r and all \tp) E 

Hi) Its convex roof extension to mixed states is given by 



E, 



it. (p) ■ = min ^2 Pi E inv (| tpi) (ipi |) , 



(1) 



where the minimum is taken over all possi- 
ble pure states ensembles that decompose p, i.e. 



lnv is an entan- 
However, we will 



The criteria above guarantee that E, 
glement monotone under LOCC 1 
explicitly show here that Ei nv is also monotone under the 
more general class of separable operations (see the upper 
bound in Theorem [!} . These measures capture genuine 
multipartite entanglement and are in general dimension- 
dependent, i.e., may not exist for some choices of the 
local dimensions di. The measures are therefore depen- 
dent on the space the states are embedded in, as we 
will see below. An example is Wootters' concurrence of 
two qubits (l5j and its higher dimension generalization 
called G-concurrence [lij]. In particular, these are the 
only SL-invariant measures for the bipartite case, with 
d\ = g?2- No SL-invariant measure exists for d\ ^= d2, 
hence a simple embedding of a 2 x 2 state into a 2 x 3 
state is not measurable. For multi-partite systems there 
is more than one measure, and their construction is based 
on SL-invariant polynomials [l4l . [l7| . An SL-invariant 
polynomial, f(\tp)(ip\), is a polynomial in the compo- 
nents of the state \ip), invariant under the action of the 
group G, i.e. f{g\^){^W) = f(U>)W), for all g E G. 
Then E inv {\ij}){^}\) = \fk{\ip) (ip\)\ 1/k is an example of an 
SL-invariant measure, where fk is a homogeneous SL- 
invariant polynomial of degree k. In particular for 3 
qubits the square root of the 3-tangle [l8| is the simplest 
multi-partite SL-invariant measure. 

In the next Lemma we prove that the convex roof 
extension on an SL-invariant measure satisfies the SL- 
invariance and homogeneity properties, similarly to the 
pure state case. This is the most important technical re- 
sult of our paper and is crucial for the proof of our main 
result, Theorem [TJ This important Lemma was over- 
looked by the community, a partial version being stated 
before in [l[ , with inequality instead of equality and valid 
only for bipartite separable operations. 

Lemma 1. Any entanglement measure defined as a con- 
vex roof extension of an SL-invariant pure state entan- 
glement measure remains SL-invariant and homogeneous 
of degree 1, i.e. 

i) E inv {gpg^) = E inv (p), for all g E G and all p E 
B{U); 

ii) Ei nv {rp) = rEi nv (p), for all non-negative r and all 
P G B{%); 



Proof. Proof of i). Let g E G and let {IV^X^il} be an 
optimal un-normalized ensemble decomposition of p, i.e. 

E inv (p)=J2 E inv(\$i)$i\)- (2) 
i 

Then {g^iji^iW} must be some ensemble decomposi- 
tion of the positive operator gpg\ in general not the op- 
timal one. Hence 

Etnvigpg* 1 ) < ^£ in „(g|^ i )('0 i |g t ) 

i 

= J2 E ^v(\A)(H) = E mv (p). (3) 



Let now {|<^fc)(</>fc|} be an optimal un-normalized ensem- 
ble decomposition of gpg^ . i.e. 



Einv{gpg^) = y^,E inv (\<t> k ){4>k\ 



(4) 



Note that {g~ l \4>k) (4>k\g~ 1 ^ } is some ensemble decompo- 
sition of p, again not necessarily the optimal one, where 
g^ 1 is the inverse of g (the inverse is guaranteed to exist 
since g has determinant 1). Then 

Einv(p) < y^-Em U (ff~ 1 |0fc)(0fc|g~ lt ) 
fe 

= ]T E mv (\^k)(4>k\) = E inv (gpgi). (5) 

k 

Now part i) follows from ([3]) and ([5]). 

Proof of ii). This follows at once, since if {\ipi)(tpi\} is 
an optimal ensemble decomposition of a density operator 
p, then {rl'i/y^il} will automatically be an optimal de- 
composition of the operator rp. Then the result follows 
from the homogeneity of Ei nv (-) for pure states. □ 

Definition 2. Let A be a completely positive trace 
preserving (CPTP) multi-partite separable oper- 
ation acting on B(Hi) <8 B^H^) ® ■ ■ ■ ® B(H n ), 
with operator-sum representation given by Kraus 



{K, 



K (n h 



with 



operators 

J2 m K^k£> <g> ■ ■ • <g> K^K^ = I <g> • ■ ■ <g> I. We 
define the entanglement resilience factor (ERF) of A to 
be 

F{A) := min ^| det^ 1 ))^ • • • | det(^))| 2 / d ", 

m 

(6) 

where the minimum is taken over all separable Kraus de- 
compositions {Km^ <g) Km^ (B> • • • <8> Km^}m of A. 

Lemma 2. ^(A) ^ 1, with equality only if A is a sep- 
arable random unitary channel, i.e. all Kraus operators 
are proportional to unitaries. 
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Proof. From the closure condition we have 

£ Kg** Kg* ® ■ ■ ■ ® ff&> = J ® ■ ■ • ® J. 



(7) 



Taking the determinant on both sides and applying the 
Minkowski's inequality for a sum of positive (semidcfi- 
nite) operators (see p. 47 of [TH) yields 



E 



det(^W)| 2 ^...|det(^))| 2 / d »<l, (8) 



where we have used that det(AA jf ) = |det(A)| 2 and 
dct(A® 5) = dct(A) fc dct(S) a for A and B square matri- 
ces of dimension ax a and 6x6, respectively. The result 
then follows at once from (|6|) and (|5J). □ 

The ERF was first introduced in Q for the particular 
case of one-sided operations of the form A^ 1 ) ® I® • • • ® J, 
where its computation reduced to an optimization prob- 
lem due to the unitary freedom in the Kraus representa- 
tion of a quantum operation [9|. However, in our case, 
it is not at all clear whether a separable operation ad- 
mits more than one separable representation (up to a 
relabelling of the Kraus operators), excluding the trivial 
cases where the operation is a tensor product of inde- 
pendent channels A = A^ ® • • • ® A^™'. For this latter 
case 



.F(A«®---® A'"') < J]>(A«), 



(9) 



since among all possible separable representations 
of A there are Kraus representations of the form 

{Kg ® ■ ■■ ® K^} ju ... dn , with £, fc Kf K W = I for 
any individual party k, and this implies that the min- 
imization (|()]) can split into n individual parts, which 
proves ©. 

The following Theorem constitutes our central result. 

Theorem 1 (Entanglement evolution). Let 
A be a multi-partite separable operation with 
operator-sum representation given by Kraus 



operators 
H m K^K$ 



{K 



(i) 



K 



(2) 



Km^}m, With 

■ ® I, let 

Ei nv (-) be an SL-invariant entanglement measure, 
and let p G B(H) be an entangled input state with 
Einv(p) 7^ 0. Then the average output entanglement is 
independent of the input state and of the entanglement 
measure, and is given by 

^7Jm PmEinv \0~m) 



Einv (p) 

]T | det(Kg)\ 2 ^ • • • | det(i4"))| 2 / d " < 1, (10) 

m 

with equality if and only if A is a sepa- 
rable random unitary channel FiZl /. Here 
Pm a rn = (Kg ® ■ ■ • ® K£>)p{K® ® ■ ■ ■ ® 

and Y^mPm = 1- 



In particular note that the upper bound in (|10[) explic- 
itly shows that E inv is an entanglement monotone under 
separable operations, as we claimed in the remarks fol- 
lowing Definition [2] 

Proof. We start by noting that 

(Kg ® • • • ® ® • • • ® i^) 1 = Pm(7m . (11) 

For the Kraus operators that have non-vanishing deter- 
minant (i.e. are invertible), we have 

|det(^))| 2 ^...|det(^))| 2 / d " 




(detKg) 1 /^ ■■■(dctKg ) ) 1 / d 



= P m o- r . 



(12) 



Note that the quantity in square brackets has determi- 
nant of magnitude one. Use now the homogeneity and 
SL-invariance property of any convex roof SL-invariant 
measure (see Lemma [T]) and apply E inv to ([T2"|) to get 

| det(Kg)\ 2 ^ ■ • • | dct(Kg^)\ 2 ' d "E mv (p) = p m E mv (a m ). 

(13) 

The Kraus operators that have determinant zero pro- 
duce states with p m Ei nv (a m ) = 0, see Lemma 1 of Q, 
and hence (fT3|) is generally valid for all Kraus operators. 
Now summing (|13[) over m yields the desired equality. 
The upper bound of 1 follows again from Minkowski's 
inequality for a sum of positive semidefinite operators 
(see the proof of Lemma [2]) . □ 

Operationally Theorem[T]can be thought of in terms of 
a separable measurement repeated multiple times, each 
time the outcome m and the resulting state a m being 
recorded and at the end the output entanglement be- 
ing averaged over all outcomes. Theorem [T] implies that 
the ratio between the average output entanglement and 
the initial entanglement can be computed explicitly as a 
function of the operator-sum representation of the sep- 
arable operation and has the following properties: i) it 
is independent of the input state; and ii) it is indepen- 
dent of the entanglement measure, as long as the latter 
is SL-invariant. 

In Theorem [T] we considered only input states with 
non- vanishing entanglement, since whenever Ei nv (p) = 
the average output entanglement must also be zero, hence 
such non-entangled states p are not of interest. 

The following Corollary then follows. 

Corollary 1. The ERF is bounded by 

E mv (A(p)) < j-j^j < Yj m P™E inv (o- m ) ^ 



Einv (p) 



Ei nv (p) 



Proof. The inequality T(A) < (J2 m PmE mv (cr m ) ) / E lnv (p) 
follows from the definition of the ERF as a minimum 
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over all Kraus decompositions, sec (j6]), and the equality 
(fTB")) proved above. 

Finally, the inequality E inv (A(p)) / E inv (p) < .F(A) can 
be proved by noting that Ei nv (h{p)) / Ei nv (p) is indepen- 
dent of the Kraus representation of the channel. We can 
therefore choose the Kraus representation that achieves 

the minimum in (JB]), {Km^ 
yields 



K 



(2) 



K, 



(n) 



} m , which 



^(A) = ^|det(^)| 



(l)\|2/di 



|det(/C j )l 



(n)\|2/d„ 



(15) 



But 

Ei 



deb{K, 



(n)\|2/d„ 



Einv (p) 

(16) 

as we just showed. Now, we use the convexity of Ei nv 
(we remind the reader that the measure Ei nv is defined 
on mixed states via a convex roof extension, see part iii) 
of Definition [1| , to get 



Ea, 



• (HP)) < EraPm £ 



Ein V (p) 



Einv(p) 



and this concludes the proof. 



(17) 



□ 



In particular, for one-sided operations, .F(A) can 
be computed explicitly [|j and, additionally, for 
one-sided operations and pure state inputs, the 
first inequality of (| 14[) is actually an equality, 
E inv (A(\il>)(iP\))/E inv (\ij)(iP\) = JT(A) 0. However, in 
the most general case of multi-partite separable opera- 
tions, this is no longer true. Consider for example a bi- 
partite separable operation A with an operator-sum rep- 
resentation given by {y/pl <8> I, ^Jl — pX g) X} } where 
< p < 1 and X = |0) (1| + |1) (0| is the usual bit- 
flip Pauli operator. In this case it is not hard to see 
that there are no other separable representations of A 
(any two Kraus representation are related by a unitary 
matrix, and it follows at once that a different repre- 
sentation of A will not be separable). The ERF J- [A] 
can then be computed explicitly, since there is no min- 
imization to be done, and equals 1. The maximally en- 
tangled state \ip) = (|00) + 1 11))/ V2 is invariant un- 
der A, hence E mv {K{\^){^\)) / E mv {\^){^\) = F[K] = 1 
and the bound is saturated. But in general an arbi- 
trary bipartite state |</>) is not invariant under A and 
has E inv (A(\(j))((j)\))/E inv (\(j))((f>\) < 1, as one can easily 
check using the closed form expression for the bipartite 
concurrence [l5j . 

Entanglement breaking. Interestingly, the capability of 
a channel to destroy entanglement is another quantity 
that is also independent of the initial state, as we will 
explain below. This is closely related to the work of [2l[ , 
in which the authors define correlation measures (classi- 
cal and quantum) based on how much "noise" one has to 
"inject" by local operations into an entangled quantum 
state to make it separable, in the limit of many copies 



and asymptotically vanishing errors. We now point out 
that such a measure can not be defined in the zero-error 
single copy case, since the quantity of interest is state- 
independent. The capability of a channel to destroy en- 
tanglement is quantified as follows. 

Definition 3. A CPTP map $ : B(Hd}_ — > B(H d ) is 
called r-partially entanglement breaking \2n] if 

SN [($ <g> I)p] ^ r, Vpe B{H d <8>H d ). (18) 

It generalizes the usual notion of entanglement break- 
ing channels [23|, the latter having r = 1. Here SN(a) 
denotes the Schmidt number 24j of the bipartite density 
operator cr, defined as 



SN(p) 



mm 

{Pi,|-0.}} 



{maxSiiflVi))}, (19) 



where the minimum is taken over all possible pure state 
ensemble decompositions of p and SR(\^pi)) denotes the 
Schmidt rank on \il>i), i.e. the number of non-zero 
Schmidt coefficients of \ipi). A state a is separable if 
and only if SN(a) = 1. 

The following proposition was recently proved in (25j 
(see Proposition 8). 

Proposition 1. A CPTP map $ : B(U d ) — > B(H d ) is 
r-partially entanglement breaking if and only if there exist 
a full Schmidt rank pure state \ip) £ T-L d ®T-L d (i.e. having 
all Schmidt coefficients strictly positive) for which 



5JV[($(8 I)\ip)(ip\] =r. 



(20) 



In other words, if a channel breaks entanglement for 
one (maximal Schmidt number) state, then it must break 
entanglement for all states, that is, its "disentangling 
power" is an intrinsic property of the channel indepen- 
dent of the input state, as long as the input has max- 
imal Schmidt number (this corresponds to the case of 
non-zero entanglement as quantified by any SL-invariant 
measure). This means that one can not use the approach 
of [2l| to quantify entanglement (or correlations) as the 
total amount of local noise one has to inject in a state to 
make it separable, since in a single copy zero-error regime 
this quantity is state independent. 

Conclusions and open guestions. We studied the evo- 
lution of mixed state entanglement under local decoher- 
ence. Our main result is summarized by Theorem [I] 
in which we proved that multi-partite mixed-state en- 
tanglement evolution under separable operations is mea- 
sure and state independent and is solely a function of 
the channel. This significantly extends the results in 
[H-Q applicable only for the restricted case of bipartite 
separable operations, and also generalize the previous 
work of [H, where only one-sided operations of the form 
AW®/ &>■••<£>/ were considered. For the first time we 
derived a closed form expression for the average entan- 
glement evolution under the most general class of separa- 
ble channels, and our results are automatically valid for 
LOCC since the later is a proper subset of the former. 
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An interesting fact is that our equality in Theorem Q] 
depends on the Kraus representation of the separable op- 
eration. Given such a separable operation, are there more 
than one separable representations of it (aside from a per- 
mutation of the Kraus operators)? A complete solution 
to this problem will allow an explicit evaluation of the 
ERF J- ( A) and provide a non-trivial upper bound on the 
entanglement evolution of arbitrary mixed states, Corol- 
lary Q] We observed numerically using non-linear opti- 
mization methods that there are many cases of separa- 
ble operations that do not admit more than one nontriv- 
ial Kraus representation, but this is far from a rigorous 
proof. 

Our results are valid for SL-invariant entanglement 
measures based on convex roof extension and are not ap- 
plicable in general to arbitrary entanglement measures. 
However, it follows at once that the convex roof extension 
is an upper bound on any other convex measure that co- 
incides with the former on pure states, i.e. both give the 
same result on a pure state. In this case Theorem [1] pro- 
vides a useful upper bound on the average entanglement 
evolution, i.e. the equality in (JTUJ) becomes an inequality. 



We discussed the "disentangling power" of a quantum 
channel and pointed out that it exhibits a similar state- 
invariance as our average entanglement ratio in Theo- 
rem [TJ This forbids a straightforward extension of [2l| to 
the single copy zero-error regime. An interesting gener- 
alization will be to define a disentangling power of a sep- 
arable operation, instead on a single-sided channel, and 
characterize the class of r-partially entanglement break- 
ing such operations. This approach may prove useful in 
extending [2l| to the single copy zero-error regime, since 
one does not expect the disentangling power to be state- 
invariant anymore. 
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